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Abstract — Physical-layer network coding improves the 
throughput of the two-way relay channel by allowing multiple 
source terminals to transmit simultaneously to the relay. 
However, it is generally not feasible to align the phases of 
the multiple received signals at the relay, which motivates the 
exploration of noncoherent solutions. In this paper, turbo-coded 
orthogonal multi-tone frequency-shift keying (FSK) is considered 
for the two-way relay channel. In contrast with analog network 
coding, the system considered is an instance of digital network 
coding; i.e., the relay decodes the network codeword and 
forwards a re-encoded version. Crucial to noncoherent digital 
network coding is the implementation of the relay receiver, 
which is the primary focus of the paper. The relay receiver 
derived in this paper supports any modulation order that 
is a power of two, and features the iterative feedback of a 
priori information from the turbo channel decoder to the 
demodulator; i.e., it uses bit interleaved coded modulation with 
iterative decoding (BICM-ID). The performance of the receiver 
is investigated in Rayeligh fading channels through error-rate 
simulations and a capacity analysis. Results show that the 
BICM-ID receiver improves energy efficiency by 0.5-0.9 dB 
compared to a non-iterative receiver implementation. 

I. Introduction 

In a two-way relay channel (TWRC), a pair of terminal 
nodes exchange information using a relay node. Each terminal 
is both a source and a destination. The terminals are assumed 
to have no direct radio link. Physical-layer network coding 
(PNC) is a technique that may be applied in the TWRC 
to increase throughput over conventional techniques by elimi- 
nating the number of transmissions necessary for information 
exchange between the terminals. The exchange is broken into 
two phases: the multiple-access (MA) phase and the broadcast 
phase. In the MA phase, the terminals transmit simultaneously 
to the relay, which receives the electromagnetic sum of the 
transmitted signals. In the broadcast phase, the relay broad- 
casts the sum of information to the terminals, each of which 
detects the other terminal's information by subtracting its own 
signal from the received sum. 

PNC may be implemented using one of two schemes: ana- 
log network coding (ANC) and digital network coding (DNC). 
In ANC, the relay forwards the received signal sum directly 
and all of the processing is performed at the terminals. While 
the benefit of ANC is a simple relay implementation, the 
disadvantage is that the noise at the relay is also forwarded to 
the terminals and the processing requirements at the terminals 
can be burdensome. In DNC, the relay performs detection of 
the network-coded bits, essentially cleaning up the noise at 
the relay. It then remodulates the signal and broadcasts to 
the terminals. The benefit of DNC is that the noise received 
at the relay is not retransmitted and the terminal receivers 
are simplified, but the disadvantage is that a more complex 
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receiver is required at the relay. Thus, a crucial aspect of 
implementing PNC is the formulation of an efficient relay 
receiver, and the selection of coded-modulation formats that 
work well with the DNC system. 

In |2|, we found that turbo-coded binary frequency-shift 
keyed (FSK) modulation is a good candidate for DNC, and 
formulated a noncoherent receiver Noncoherent FSK is useful 
in scenarios where the received signal phase is corrupted and 
cannot be tracked accurately. Examples of phase corrupting 
effects include imperfect or unsynchronized oscillators and 
Doppler shift. Tracking phase is even more challenging for 
a DNC relay receiver than for a conventional point-to-point 
receiver, as the relay must track two phases simultaneously 
rather than one. Noncoherently detecting the FSK signal elim- 
inates the phase tracking requirement. Multi-tone orthogonal 
FSK modulation is beneficial when energy efficiency has 
higher priority than bandwidth efficiency. If £b/No is held 
constant, then increasing the modulation order increases the 
energy per symbol, which increases the minimum distance 
between symbols, leading to lower error rates. However, 
bandwidth usage is proportional to the modulation order, 
implying a trade-off between bandwidth and energy efficiency. 

A particular technique for combining binary channel coding 
and A/-ary modulation with i\f > 2 is bit-interleaved coded 
modulation (BICM) |3|. A binary channel codeword is gen- 
erated, interleaved, and passed to an A/-ary modulator, which 
maps codeword bits to symbols for channel transmission. The 
receiver demodulates the symbols, producing soft estimates 
of each bit. A binary soft-input channel decoding scheme 
is applied to the soft estimates. The energy efficiency of 
BICM can be improved by feeding back information from 
the channel decoder to the demodulator, and performing 
iterative detection. The notion of feeding back information 
from decoder to demodulator is called BICM with iterative 
decoding (BICM-ID) [4J. 

The present work extends p) by considering turbo-coded 
multi-tone FSK modulation, and deriving an iterative receiver 
capable of performing iterative AI — ary demodulation and 
turbo-decoding. The key contribution of this work is a DNC 
relay receiver capable of noncoherent operation which utilizes 
BICM-ID for improved performance. BICM-ID improves the 
error rate performance of the relay receiver versus BICM, 
while adding computational complexity to the receiver, intro- 
ducing a trade-off in receiver design. A receiver formulation 
considering BICM-ID and FSK modulation for single-source, 
point-to-point channels is given by 151. The present work 
extends the receiver in |5| to support the two-source TWRC 
model and digital network coding. 

A Turbo BICM channel decoding scheme for the relay in 
the MA phase of DNC is given by |6|. In contrast to the 
current work, the scheme in |6| considers only coherent recep- 
tion and two-dimensional modulation. Further, this work does 
not consider the BICM-ID feedback scheme from decoder to 
demodulator |7| gives a channel model for the DNC MA 
phase that considers symbol and frame asynchrony between 



transmissions by the terminals. Techniques for optimal de- 
tection of network-coded information at the relay employing 
LDPC channel coding and BPSK modulation under symbol 
and frame misalignment are presented. To our knowledge, no 
prior work has considered the application of BICM-ID to a 
noncoherent DNC relay receiver. 

In general, there are several ways to implement channel- 
coded DNC depending on the specific order in which de- 
modulation, channel and network coding are applied at the 
relay. In |8|, the performance of different channel-coded DNC 
techniques are compared. The technique in the present work, 
in which the relay demodulator forms likelihood ratios of the 
network-coded bits for channel decoding, is contrasted with a 
technique in which the demodulator performs decoding using 
probability mass functions of the arithmetic sum of received 
symbols at the relay. It is shown that the technique utilizing 
arithmetic sums achieves superior capacity to bit-level network 
coding scheme in particular SNR regions, however, the spe- 
cific channel code applied to arithmetic network coding is a 
repeat-accumulate code designed specifically for the scheme, 
and the extension to codes such as LDPC or Turbo requires 
modifications to existing channel coding algorithms, which 
have not yet been considered. Our usage of the network-coded 
bit scheme is based on the assumption that its implementation 
in existing systems requires no modification to existing chan- 
nel coding schemes, which is a practical advantage. 

An outline of the remainder of this paper is as follows. Sec- 
tion |ll] develops the system model used throughout the paper. 
Section III presents the development of the relay demodulator 
capable or performing digital-network coded BICM-ID at the 
relay in the TWRC. Section IV presents simulated error rate 
performance and capacity analysis of the relay receiver. 

II. System Model 

This section presents the system model used throughout the 
work. Modulation, channel coding, channel model, and relay 
reception are described in precise detail. An overview of the 
iterative decoding process at the relay is given. The system 
model illustrating transmission by the terminal nodes to the 
relay is shown in Fig. [T| 

A. Transmission by Terminal Nodes 

The terminal nodes A/i, i € {1,2} generate binary in- 
formation sequences — [wi i, i] having length K. 
A rate-rs turbo code is applied to each u^, generating a 
length L = Kjrg binary channel codeword, denoted by b'^ = 
[5i i, 6L,i]- The codeword is passed through an interleaver, 
modeled as a permutation matrix 11 having dimensionality 
i X i: b., = b^n. Let V = {0, ...,M - 1} denote the set 
of integer indices corresponding to each FSK tone, where M 
is the modulation order. The number of bits per symbol is 
/i — log2 M. The codewords b^ at each node are divided into 
Nq ~ L/fi sets of bits, which are passed to an Af-ary FSK 
modulator. The modulator maps each set to an M -ary symbol 
qk,i € 25, where k denotes the symbol number, and i denotes 
the terminal. The modulated signal transmitted by terminal Mi 
during signaling interval kTg < t < {k + l)Ts is 
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where Sfe i(i) is the transmitted signal, fi is the carrier 
frequency of terminal JVi, and Tg is the symbol period. 

The continuous-time signals Sk,i{t) are represented in dis- 
crete time by the set of column vectors {e^j. . : qk^i € V}. 
The column vector e^^. . is length M, contains a 1 at vector 
position Qk.i, and elsewhere. The modulated codeword 



from terminal A/i is represented by the matrix of symbols 
Xi [xi.i, XAr_j.i], having dimensionality AI x Nq, where 

B. Channel Model 

All channels are modeled as flat-fading channels having 
independent gains for every symbol period. The complex- 
valued channel gain from node Mi to the relay during a 
particular signaling interval k is denoted by hk.i- The gain 
is represented as h^^i = ak,ie'^^-\ where ak,i is the received 
amplitude and O^.i is the phase, which depends on the phase 
shift of the channel and corruption in the reference signals 
at the terminal and relay induced by hardware imperfections. 
The amplitudes of the gains are selected such that the received 
energy at the relay from node Mi is 8i 



(2) 



Consider transmission of a single frame of Nq symbols to 
the relay. The received frame is 

Y = XiHi + X2H2 + N (3) 

where is an Nq x Nq diagonal matrix of channel co- 
efficients having value h^^i at matrix entry (n, n) and 
elsewhere, and N is an M x A'^ noise matrix. A single pair 
of channel-corrupted symbols received at the relay is denoted 
by y, and defined as a channel observation. In terms of this 
definition, Y = [yi, yAr^^], where y^ denotes the fc-th 
channel observation. Denote the fc-th column of N by n^. 
Each column is composed of zero-mean circularly symmetric 
complex Gaussian random variables having covariance matrix 
A'qIm; i-e., Hfc ^ Mc{0, NqIm)- is the one-sided noise 
spectral density, and Im is the M-by-M identity matrix. 

C. Relay Reception 

The goal of the relay receiver is to detect the network-coded 
combination of information bits transmitted by the terminals, 
u — Ui U2. The relay receiver takes as input the frame 
of channel observations Y. The symbols transmitted by the 
terminals are assumed perfectly synchronized at the receiver. 
Demodulation and iterative channel-decoding are applied to 
the received frame to detect u. Define the network codeword 
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Fig. 1. System Model - TWRC DNC MAC Phase 



Since the turbo code is a systematic linear code, b forms a 
code from the codebooks used by the terminal nodes, thus, 
the channel decoding operation yields a hard decision on the 
network coded message bits u. 

The demodulator and decoder are implemented using the 
soft-input soft-output algorithm described in |9|. The de- 
modulator takes as input the matrix of received symbols Y 
and a-priori probability (APP) information v and produces 
extrinsic information z, v = [vi,...,vl], z = [zi,.-.,zl]. 
Precise description of the quantities and is provided 
in Section IIIll The extrinsic information is deinterleaved to 
produce z' = zU^ and passed to the decoder The decoder 
refines the estimate of z', producing new extrinsic information 
v' which is interleaved to produce v = v'll and returned to 
the demodulator 

The average symbol signal-to-noise ratio transmitted by 
each terminal Si/No i E {1,2} is known at the demodu- 
lator. The demodulator may operate using several cases of 
knowledge of the channel gains hk^i (termed channel state in- 
formation or CSI): the case in which the gains are completely 
known, the case in which only the fading amplitudes ak.i are 
known, and the case in which no information about the gain 
is known (no CSI). 

Specific details of the demodulator are discussed in Section 
Hm The details of the channel decoder have been treated 
extensively in the literature |9|, 1 10|, and will not be discussed 
here. Note that a feed-forward BICM receiver without BICM- 
ID does not feed back extrinsic information from the decoder 
to the demodulator 

The relay encodes and modulates the network-decoded 
message bits u and broadcasts its modulated signal to the 
terminal nodes. The signal traverses two independent fading 
channels, and the terminal nodes receive independently faded 
versions of the message bits corrupted by white noise. The 
terminal nodes demodulate and decode the signal received 
from the relay to form estimates of u. Let u denote the 
detected sequence at Afi and u denote the detected sequence at 
JV2- Each terminal estimates the information bits transmitted 
by the opposite node by subtracting its own information 
sequence from the sequence detected from transmission by the 
relay: g2 = u Ui at Afi and gi = u ® U2 at Af2- Since the 
links from the relay to the terminals are conventional point-to- 
point links with no interfering transmissions, specific details 
of the terminal-node receivers are omitted. 

III. Soft N-FSK Relay Demodulator 

The soft-output relay demodulator maps the received sum 
of symbols from the terminals to log-likelihood ratios of 
the network-coded bits. The demodulator operates iteratively, 
using extrinsic information fed back from the channel decoder 
to refine the information prior to each decoding iteration. After 
a specified number of iterations has been reached, the decoder 
makes a hard decision on the network-coded bits. 

The demodulator processes a frame of channel observations 
Y one observation at a time. Since the operation performed 
on each observation is the same, we may drop the dependence 
on a particular signaling interval in the frame to simphfy 
the notation. Denote a single received channel observation 
as y. During the first demodulation and decoding iteration, 
the demodulator computes the probability of every possible 
combination of symbols transmitted by the terminals: P{q; I), 
where q is defined as the tuple 

<? = 92) qi,q2eV qeVxV (5) 

and qi and (72 denote the symbols from terminal JVi and A/2, 
respectively. We will refer to (7 as a super-symbol. Define this 
operation as the super-symbol probability mapping stage. The 



symbol probabilities P{q; I) are fixed for all demodulation 
and decoding iterations. Since the cardinality of 2? x 2? is 
A/^, the relay receiver computes AP probabilities, versus a 
conventional point-to-point receiver which only computes M 
probabilities, as only one terminal is present. 

On the first and subsequent decoding iteration, the symbol 
probabilities are transformed to the set of /i log-likelihood 
ratios associated with each network-coded bit mapped to 
the super-symbol. Denote this operation as digital network- 
coded soft mapping (DNC-SOMAP). A general description 
of SOMAP for the point-to-point channel is given by |9|. The 
DNC-SOMAP takes as input the symbol probabilities F{q; I) 
and extrinsic information represented by bit probabilities as- 
sociated with each network coded bit P{c; I) fed back from 
the channel decoder, where c denotes the ji network coded 
bits mapped to q, P{c;I) = {P(cfe;/),0 < fc < ^ - 1}, 
Cfc = Ci © C2 fc, and q ^ denotes the fc-th bit mapped to 
symbol q^. i G (l,2}. On the first demodulation iteration, no 
decoding has been performed, and the bit probabilities are as- 
sumed equally likely. The DNC-SOMAP produces estimated 
probabilities of values taken by c: P{c; O) = {P{ck; /), < 
k < 1.1- 1}. 

The input distributions to the DNC-SOMAP with respect to 
the super-symbol are represented as probabilities, while the in- 
put and output distributions with respect to the network-coded 
bits mapped to each super-symbol are represented as log- 
likelihood ratios. Log-likelihood representation facilitates soft- 
decision decoding. The input representing the log-likelihood 
ratio of the rn-th bit mapped to the super-symbol is related to 
the input distribution by 

Vk=\og l\'"'^l''^\ , 0<fc<M-l. (6) 
P(cfc = 0; /) 

The output representing the log-likelihood ratio of the k- 
th bit mapped to the super-symbol is related to the output 
distribution by 

^fe-log pi''-^^'^i , 0<k<^^-l. (7) 
P(cfc = 0; O) 

The DNC-SOMAP output distribution is related to the input 
distributions by 

Substituting the specific values of the distribution (j6|l into the 
expression for output ([8]l, 

P(c,=£;0)= Y: Piy\l)li^ (9) 

q:ck=e j=o 

The output log-likelihood ratio of the DNC-SOMAP may 
be found by combining ([9]l and (|7]l: 

q:Ck = l j=0 

zk = log '-^^ (10) 

q:Cfc=0 j=0 



where the term (1 + e""^) cancels in the ratio. For the pur- 
pose of numeric implementation, it is useful to simplify this 
expression using the max-star operator 



max*{xi} = log ■ 



(11) 



where the binary max-star operator is max*(a::,y) = 
max(a;,y) + log(l + e^'^^^l) and multiple arguments im- 
ply a recursive relationship; for example: max*(a;,y, 2; 



max max z)). Applying the max-star operator to ( 10 1 
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A non-iterative BICM demodulator does not utilize feedback 



from the decoder, so it is implemented using (12 1 setting all 



0. The values taken by the pdf p{y\q) are dependent on 
the available channel state information. Description of these 
pdfs is given in the following subsections. 

A. Noncoherent Reception with CSI 

Considering noncoherent reception with CSI, the phases of 
hi and /i2 are not known. The pdf of the super-symbol p{y\q) 
takes different forms depending on the values of the symbols 
Qi and (72- When qi and q2 are different, the pdf takes the 
form |2j) 



p{y\q) = exp 
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where Jq denotes the zeroth-order modified Bessel function of 
the first kind, and j/^^ and t/^^ denote the q^ and (72-th entries 
of the received channel observation y. When qi and q2 are 
the same, the form of the pdf is [2] 



p{y\q) = exp 
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where a is approximated as a = vaf + af . A discussion of 



this approximation is found in |2 
mapper computes the values of 



The symbol probability 
and 



13 1 and (14i which are 
substituted into (fT2]i to compute the output log-ukelihood ratio 
of the DNC-SOMp. 

B. Noncoherent Reception Without CSI 

When the relay possesses no knowledge of the channel 
gains, the pdf p{y\q) is a function of the average symbol 
energy transmitted by the terminals and the noise variance. 
In the case that the symbols qi and (72 are different, the pdf 
takes the form ||TTJ 
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Fig. 2. Error-rate performance of the relay demodulator as a function 
of modulation order and decoding iterations. The relay receiver possesses 
CSI. The data sequence is length K = 1229. The channel code rate is r = 
1229/2048 0.6. Solid lines denote BICM, and dashed lines denote BICM- 
ID. With each modulation order, from right to left, the number of decoding 
iterations is 1, 10, and 30. 



and when qi and (72 are the same |[TT] 
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The symbol probability mapper computes the values of ( [T5[ ) 
and (fl6| which are substituted into (fT2| to compute the output 
log-lkelihood ratio of the DNC-SOMAP. 

IV. Simulation Study 

In this section, simulated error rate performance and ca- 
pacity analysis is shown considering the demodulator derived 
in Section |lll] under the channel model described in Section 
Inj Error rate and capacity is simulated for the MAC phase 
or DNC two-way relaying only, as the the broadcast phase is 
equivalent to a pair of point-to-point links, which has been 
thoroughly discussed in the literature. 

Error rate performance is simulated as a function of FSK 
modulation order, channel state information at the relay, 
channel decoding iterations, and presence or absence of in- 
formation feedback from decoder to demodulator via BICM- 
ID, and signal-to-noise ratio. All error rate simulations utilize 
soft-decision channel coding. The specific channel code is the 
Turbo code described by the UMTS standard |12|. Channel 
capacity is simulated as a function of channel state information 
and modulation order The results of simulation are interpreted 
to provide recommendations for relay receiver configuration. 

A. Error Rate Performance 

In this subsection, the specific results of error rate simula- 
tion are presented. All simulations utilize a UMTS Turbo code 
with rate K/L = 1229/2048 w 0.6. The FSK modulation 
orders utilized at the terminals are M ~ {2, 4, 8}. The demod- 
ulator is simulated with and without ch annel s tate information, 
as described in subsections IIII-AI and IIII-BI The number of 
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Fig. 3. Error-rate performance of the relay demodulator as a function of 
modulation order and channel state information. The number of decoding 
iterations for all curves is 30. The data sequence is length K = 1229. The 
channel code rate is r = 1229/2048 ^ 0.6. Solid lines denote BICM, and 
dashed lines denote BICM-ID. 
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Fig. 4. Error rate performance of the relay demodulator as a function of 
decoding iterations and channel state information. The modulation order is 
M = 4. The data sequence is length K = 1229. The channel code rate is 
r = 1229/2048 0.6. Solid lines denote BICM, and dashed lines denote 
BICM-ID. Within each case of channel state information, from right to left, 
the number of decoding iterations is 1,2,4,30. 



decoding iterations is chosen from {1, 2, 4, 30}. The range of 
decoding iterations was chosen between 1 — 30 as negligible 
performance improvement was observed beyond this range 
for the selected system and channel parameters. Decoding is 
performed with BICM and BICM-ID for modulation orders 
M > 2. 

Error rate performance of the relay demodulator as a 
function of modulation order and decoding iterations is shown 
in Fig. [2] Comparing curves utilizing the same number of 
decoding iterations, increasing M from 2 to 4 and 4 to 
8 improves energy efficiency by 4 and 2 dB respectively, 
regardless of the number of decoding iterations or decoder 
feedback. The performance of the cases of M = 2 with 
BICM utilizing both 10 and 30 decoding iterations has better 
energy efficiency than the case of AI = 4, BICM, and 1 
decoding iteration. Similar behavior is observed in comparing 
M = 4 to i\/ = 8. This result implies an energy efficiency 
trade-off - particular error rate and energy operating points 
may be achieved by varying AI or the number of decoding 
iterations. The absolute performance improvement of BICM- 
ID vs. BICM for AI = 4 and 8 is approximately 0.5 
and 0.7 dB, respectively, implying that relative performance 
improvement increases with modulation order. 

Fig. |3] shows error rate performance of the relay demod- 
ulator as a function of modulation order and channel state 
information. The number of decoding iterations is fixed at 
30. The relative performance as a function of CSI is roughly 
the same between curves having the same modulation order 
and decoder feedback. The absolute performance improvement 
of BICM-ID over BICM considering AI = 4 is about 0.75 
and 0.5 dB for no CSI, and CSI, respectively. For AI = 8, 
the performance improvement is about 0.9 and 0.6 dB. The 
implication is that BICM-ID provides a greater relative im- 
provement in the absence of channel state information and 
with increasing AI. A notable result is that BICM-ID with no 
CSI has better energy efficiency than BICM with CSI. This 
implies that decoder complexity may be closely traded off 
with channel knowledge to achieve a given energy efficiency. 

Error rate performance as a function of decoding iterations 



and channel state information is shown in Fig.j4] The modu- 
lation order is fixed at AI = 4. Considering BICM curves, the 
performance difference between CSI and no CSI for a given 
iteration decreases as the number of iterations is increased. 
In the case of 1 iteration, the performance difference is about 
2 dB, while for 30 iterations, the difference is about 0.5 dB. 
The same property holds for BICM-ID. BICM-ID allows the 
no CSI cases to outperform the CSI cases by varying the 
number of decoding iterations. Consider the case of no CSI, 
BICM-ID and 4 decoding iterations. This case outperforms the 
CSI cases using both BICM and BICM-ID using 2 iterations. 
This example clearly illustrates the potential design trade-off 
in utilizing decoder iterations and channel state information. 
Likewise, for 30 iterations, the no CSI case with BICM-ID 
outperforms all CSI cases except the case also utilizing 30 
iterations and BICM-ID. 

B. Binary Information Rate 

This subsection presents simulated values of the binary 
information rate of the relay demodulator as a function of 
modulation order and channel state information. Binary infor- 
mation rate is a measure of maximal throughput achievable for 
a particular BICM receiver configuration and channel model 
as a function of signal to noise ratio. The metric can also 
be interpreted as the minimum energy required for error-free 
when a capacity approaching code is used. 

Binary information rate is computed by ||5| 



R = 1 



log2(e) 



^£;[max*(0,A(6fe)(-l)'"=] (17) 



fc=0 



where is the fc-th network-coded bit associated with a 
particular super-symbol, and A(fefc) is the lo g-lik elihood ratio 
associated with the bit. The expectation in ( 17 1 is computed 
by Monte Carlo simulation. 

A sequence of bits having length K = 10000 are generated 
by each user, divided into groups having length /i. Each group 
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Fig. 5. BICM capacity at output of relay demodulator as a function of 
modulation order and channel state information. Dashed and solid lines denote 
no CSI and CSI, respectively. For reference, the Ei,/N(i value required to 
reach an eiTor rate of 10"'' for several feedback configurations and channel 
states are shown. All receivers perform 30 decoding iterations. 



is mapped to symbol, modulated, and transmitted over the 
channel described in Section |ll] The bitwise LLRs of the 
corresponding super-symbols are computed at the relay using 
the demodulator described in Section (|III]l. The LLRs and 
network coded bits in the frame are substituted into ^TT) to 
compute the value of R. The expectation in ( [T7| ) is computed 
over one million trials and for several values oFsymbol energy 
to noise ratio Es/Nq. The bit energy to noise ratio can then 
be computed as £i, /Nq — £s/Nq/R. 

Simulated values of binary information rate are shown in 
Fig. [3] From top to bottom, the dashed lines and solid lines 
represent binary information rate with no CSI and CSI respec- 
tively, for modulation orders M ~ {2, 4, 8}. Energy efficiency 
improves with increasing modulation order, as expected using 
orthogonal modulation. The gap between the energy efficiency 
with and without CSI increases with modulation order. A 
design implication is that channel estimation may be more 
beneficial in systems utilizing higher modulation orders. The 
minimum of each curve represents the most energy-efficient 
rate for a given level of CSI and modulation order. For M = 
2, 4, 8, the most energy efficient rates are approximately 
0.2, 0.26, and 0.3, respectively. The design consideration is 
that higher throughputs are achievable using less energy as 
modulation order is increased. 

The rate curves shown in Fig. [5] represent the achievable 
energy efficiency for a given code rate. Actual systems will 
operate at efficiencies which are suboptimal. Improvements 
in system design yield improvements in energy efficiency. 
In this work, the application of feedback from demodulator 
to decoder using BICM-ID improves efficiency over BICM 
which does not utilize feedback. 

Several points representing the Eb/No required to reach 
a simulated error rate of 10^^ are shown on Fig. [s] For a 
particular rate and modulation order, points are shown for 
simulation cases considering different receiver configurations 
and CSI levels. For all modulation orders, we see that the effi- 



ciency difference between CSI and no CSI receivers increases 
with code rate. This implies that channel estimation yields 
higher performance improvement as code rate increases. For 
M = {4, 8}, BICM-ID improves efficiency over BICM, with 
higher performance improvement as code rate increases. 

Considering the trade-off between decoder feedback and 
channel state information at a particular modulation order, 
we observe that at code rates B = 0.4 and R — 0.6, 
a system applying BICM-ID with no CSI outperforms a 
system applying BICM with CSI. This implies that decoder 
feedback may be traded off with channel state information 
as a means of performance improvement at low rates. At rate 
R = 0.9, considering the configuration of BICM with no CSI, 
applying channel state information yields more performance 
improvement than BICM-ID. 

V. Conclusion 

This work presents a relay receiver capable of performing 
digital network coding in the two-way relay channel using 
soft-input channel decoding and feedback from decoder to 
demodulator, termed BICM-ID. Simulation results using the 
UMTS Turbo code, 4, and 8-ary modulation, and different lev- 
els of channel state information show error rate improvements 
between 0.4-0.9 dB over non-BICM-ID systems. A natural 
extension of this work is adapting the relay receiver to support 
coding schemes that perform channel decoding on the sums 
of received symbols, rather than network-coded bits, which 
has been shown to improve capacity |8j. 
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